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Abstract 

We present a detailed spectroscopic study of the hot gas in the Galactic halo toward the direction of a 
blazer PKS 2155-304 (z =0.117). The O VII and O VIII absorption lines are measured with the Low and High 
Energy Transmission Grating Spectrographs aboard Chandra, and the OVII, OVIII, and NelX emission 
lines produced in the adjacent field of the PKS 2155-304 direction are observed with the X-ray Imaging 
Spectrometer aboard Suzaku. Assuming vertically exponential distributions of the gas temperature and 
the density, we perform a combined analysis of the absorption and emission data. The gas temperature 
and density at the Galactic plane are determined to be 2.5(+0.6, -0.3) x 10 6 K and 1.4(+0.5, -0.4) x 10~ 3 
cm~ 3 and the scale heights of the gas temperature and density are 5.6(+7.4, — 4.2) kpc and 2.3(+0.9, — 0.8) 
kpc, respectively. These values are consistent with those obtained in the LMC X-3 direction. 
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1. Introduction 

X-ray observations of edge-on spiral galaxies revealed 
the existence of hot gas at temperatures of ~ 10 6 K ex- 
tending a few kpc beyond the disk (e.g. Wang et al. 
2001; Wang et al. 2003; Strickland et al. 2004; Li et al. 
2008; Yamasaki et al. 2009). The origin of energy and ma- 
terial in such a hot halo has not been clarified. Feedback 
from supernovae (SNe) as galactic wind or fountain and 
heated primordial gas are possible candidates (Norman & 
Ikcuchi 1989). In any cases, halo gas plays important roles 
in galactic evolution through chemical circulation and in- 
teraction between galaxies and the intergalactic medium. 

The hot gaseous halo in and around the Milky- Way has 
been investigated for a long time. For instance, ROSAT 
All Sky Survey (RASS) quantitatively mapped the spatial 
distribution of the Soft X-ray Background emission (SXB; 
Snowdcn et al. 1997). The Cosmic X-ray Background 
(CXB) component extrapolated from the discrete hard 
X-ray sources could explain only about half of the SXB, 
leaving the soft X-ray emission below 1 keV being of a dif- 
fuse origin. With the high resolution X-ray microcalorime- 
ter flying on a sounding rocket, McCammon et al. (2002) 
detected emission lines of hydrogen- and helium-like oxy- 
gen, neon, and iron ions from about 1 steradian of the 
sky, which suggests that the emitting gas is of a thermal 
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nature and at temperatures of T~ 10 6 K. The existence 
of the hot gas in and around the Milky- Way is consistent 
with the Chandra observations of nearby edge-on spiral 
galaxies. However, because these emission data carry very 
little distance information, the properties of the global hot 
gas, like its density, temperature, and their distributions, 
are still poorly understood. 

A combined analysis of high resolution absorption and 
emission data provides us with a powerful diagnostic of 
properties of the absorbing/emitting plasma. Absorption 
lines measure the column density of the absorbing mate- 
rial, which is an integration of the density of the absorbing 
ions along a sight line. In contrast, emission line inten- 
sity is sensitive to the emission measure, which is propor- 
tional to the density square of the emitting plasma. Thus, 
a combination of the emission and absorption data natu- 
rally yields the density and the size of the corresponding 
absorbing/emitting gas. 

With significantly improved spectral resolution of cur- 
rent X-ray instruments, we are now able to observe the 
needed high resolution absorption and emission lines pro- 
duced in the hot plasma. For instance, the X-ray absorp- 
tion lines at z = 0, in particular the helium- and hydrogen- 
like O VII and O VIII lines, are detected in spectra of many 
galactic and extragalactic sources (e.g. Futamoto et al. 
2004; Yao & Wang 2005; Williams et al. 2007). Recently, 
Fang et al. (2006) and Bregman & Lloyd-Davies (2007) 
find that the O VII absorption line can always be de- 
tected in an AGN spectrum as long as the spectrum is of 
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high signal-to-noise ratio. On the other hand, the X-ray 
Imaging Spectrometer (XIS) aboard Suzaku can also re- 
solve emission lines produced in a diffuse emitting plasma 
at temperatures of T~ 10 6 K. And indeed, the O VII and 
OVIII lines have been detected in nearly all directions 
(e.g., Smith et al. 2007; Shelton et al. 2007). Recently, 
a systematical study of emission lines of the hot gas in 
and around the Galaxy has been conducted by Yoshino 
et al. (2009), who report the OVII and OVIII lines in 14 
blank sky observations with the XIS and conclude that 
the linc-of-sight mean temperatures of the emitting gas 
has a narrow distribution around 2.3 x 10 6 K. Since the 
ion fractions of O VII and O VIII and their K-transition 
emissivities are very sensitive to gas temperature at ~ 10 6 
K , a combined analysis of these emission and absorp- 
tion lines will also constrain the gas temperature and its 
distribution without the complexity of relative chemical 
abundances of metal elements. 

Although this combined analysis method has long been 
applied in the ultraviolet wavelength band (Shull & Slavin 
1994), its application in the X-ray band just began. 
Complementing the high resolution absorption data ob- 
served with Chandra with the broadband emission data 
obtained with RASS, Yao & Wang (2007) firstly at- 
tempted to conduct the combined analysis in the X-ray 
band to infer the hot gas properties in our Galaxy. They 
also proposed a model for the Galactic disk assuming the 
temperature and density of the hot gas fading off expo- 
nentially along the vertical direction. They concluded that 
the OVII and OVIII absorption lines observed along the 
Mrk 421 sight line are consistent with the Galactic disk 
origin. Yao et al. 2009 further constrained this disk model 
by jointly analyzing the high resolution absorption data 
obtained with Chandra along the LMC X-3 sight line and 
emission data observed with Suzaku in the vicinity of the 
sight line. They estimated gas temperature and density 
at the Galactic plane and their scale heights as 3.6 (+0.8, 
-0.7) xlO 6 K and 1.4 (+2.0, -1.0) xl0~ 3 cm~ 3 and 1.4 
(+3.8, -1.2) kpc and 2.8 (+3.6, -1.8) kpc, respectively. 
These results are consistent with the early findings by Yao 
& Wang (2007), i.e., the SXB can be explained by a kpc- 
scale halo around our Galaxy. 

In this paper, we present the second case study of the 
combined analysis of high resolution absorption and emis- 
sion lines. The absorption lines are observed with Chandra 
along a blazer, PKS 2155-304 sight line and the emission 
lines are obtained with Suzaku observations of the vicinity 
of the sight line. In Section 2, we describe our observa- 
tions and data reduction process. We perform our data 
analysis in Section 3 and discuss our results in Section 4. 

2. Observations and Data Reduction 

2.1. Chandra Observations and Data Reduction 

Chandra observed PKS 2155-304 many times. There 
are two grating spectrographs (the low and high energy 
transmission grating spectrographs; LETG and HETG) 
and two sets of detectors ( the advanced CCD imag- 
ing spectrometer; ACIS and the high resolution camera; 
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Fig. 1. RASS 0.1 - 2.4 keV band X-ray map in the vicinity of 
PKS 2155-304 (the bright source at the center) and the XIS 
field of view of the two presented observations. 

HRC) aboard Chandra 1 . In this work, we used all ob- 
servations available to the date of 2009 March, except for 
some observations made with non-standard configuration 
of ACIS (i.e., putting source outside the CCD-S3 chip) 
to avoid spectral resolution degradation. The data used 
in this work include 46 observations with an accumulated 
exposure time of 1.07 Ms. 

We followed the standard scripts to calibrate the ob- 
servations 2 . When extracting the grating spectra and 
calculating the instrumental response files, we used the 
same energy grid for all observations with different grat- 
ing instruments and/or with different detectors for ease of 
the adding process described in the following. For those 
HETG observations, we only use the first order grating 
spectra of the medium energy grating (MEG) to utilize its 
large effective area at lower (< 1 keV) energy. For those 
observations taken with the HRC, we further followed the 
procedure presented in Yao et al. (2009) to extract the 
first order spectra of the LETG. We then added the first 
grating order spectra of all observations to obtain a sin- 
gle stacked spectrum and a corresponding instrumental 
response file. 

2.2. Suzaku Observations and Data Reduction 

We observed the emission of the hot diffuse gas toward 
two off-fields of the PKS 2155-304 sight line during the 
A02 program (Table 1). To minimize confusion by stray 
lights from the PKS 2155-304 and to average out the pos- 
sible spatial gradient of the diffuse emission intensity, the 

1 please refer to the Chandra Observatory Guide for more infor- 
mation: 

http: / /cxc. harvard.edu/proposer/POG/html/indcx. html 
Please refer to the CIAO script for more information: 
http: / / cxc.harvard.edu/ciao / guides / 
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Table 1. Suzaku Observation Log 

Szl Sz2 

(a, 5) in J2000 (degrees) (329.2236,-30.5193) (330.1731,-29.9560) 

(£,b) in Galactic coordinate (degrees) (17.1809,-51.8544) (18.2418,-52.6081) 

Observation ID 503082010 503083010 

Observation start times (UT) 18:32:39, 2008 Apr 29 08:31:41, 2008 May 2 

Observation end times (UT) 08:30:08, 2008 May 2 17:30:19, 2008 May 4 

Exposure time 90ks 87ks 

Exposure after the data reduction 51. Iks 56.3ks 
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Fig. 2. (a) XIS light curve in 0.3-2.0 keV and (b) solar wind 
proton flux calculated using the data of ACE SWEPAM in 
Szl (top) and Sz2 (bottom) observation periods. The time is 
plotted from the beginning of each observation with Suzaku. 
The time bin of proton data are shifted 5000 seconds to correct 
for the travel time of the solar wind from the ACE satellite 
to the Earth. The dashed lines in the bottom panel indicate 
the threshold of the proton flux as 4 X 10 s cm~ 2 s — 1 . 

two fields were chosen to be 30' away from the PKS 2155- 
304 and in nearly opposite directions (Fig. 1). With this 
configuration and the roll angle of the XIS field of view, we 
estimate that stray lights from PKS 2155-304 contribute 
no more than 10% to the observed X-ray emission in 0.3- 
1.0 keV energy range. Our observation pointings are away 
from the southern edge of Radio Loop I (Bcrkhuijsen et 
al. 1971). Thus we consider that there is no contamina- 



tion of the emission from the Loop I in our observations. 
This is supported by the observational results that there 
is no EUV enhancement in this direction (Sembach et al. 
1997). 

Our observations were taken with the CCD camera X- 
ray Imaging Spectrometer (XIS; Koyama et al. 2007) on 
board Suzaku (Mitsuda et al. 2007). The XIS was set to 
the normal clocking mode and the data format was either 
3 x 3 or 5 x 5, and the spaced-raw charge injection (SCI) 
was applied to the data during the observations. We used 
processed data version 2.2.7.18 for the two observations. 
In this work, we only used the spectra obtained with XIS1. 
Compared to the other two front side-illuminated CCDs, 
XIS0 and XIS3, XIS1 is a backside-illuminated CCD chip 
and is of high sensitivity at photon energies below 1 keV. 
We found no point sources in the FOV, thus we used the 
full CCD field of view in further analysis to increase the 
photon counts because X-ray from the calibration sources 
do not affect the soft X-ray spectrum below 5 keV. 

We adopted the standard data selection criteria to ob- 
tain the good time intervals (GTIs), i.e. excluding expo- 
sures when the line of sight of Suzaku is elevated above the 
Earth rim by less than 20° and exposures with the "cut-off 
rigidity" less than 8 GV. We checked the column density 
of the neutral oxygen in the Sun-lit atmosphere in the line 
of sight during the selected GTIs, and excluded the expo- 
sures when the column density is larger than 1.0 x 10 15 
cm~ 2 to avoid significant neutral oxygen emission from 
Earth's atmosphere (Smith et al. 2007). We created X-ray 
images in 0.4-1.0 keV energy range for the two observa- 
tions, and found no obvious discrete X-ray sources in the 
fields. 

In the last step, we excluded those events severely con- 
taminated by the X-ray emission induced by the solar- 
wind charge exchange (SWCX) from geocorona (Fujimoto 
et al. 2007), meeting either of the following two crite- 
ria by Yoshino et al. (2009). The first one is the solar 
wind flux (Fig. 2). We used the solar wind data obtained 
with the Solar Wind Electron Proton and Alpha Monitor 
(SWEPAM) aboard the Advanced Composition Explorer 
(ACE ) and removed the time intervals when the proton 
flux in the solar wind exceeds 4 x 10 s cm -2 s _1 (Masui 
et al. 2009). ACE is in LI of the Solar-Earth system, 
1.5xl0 6 km away from the Earth and assuming average 
solar wind velocity as 300 km s _1 , we corrected traveling 
time of the solar wind from LI to the Earth. The second 
criteria is the Earth-to-magnetopause (ETM) distance in 
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the line sight of Suzaku (Fujimoto et al. 2007), which is 
required to be > 5Re- We found that about 20% and 
5% of the exposure time of our 1st and 2nd observations 
meets the first criteria and no time meets the second cri- 
teria. Thus we exclude that 20% and 5% from 1st and 
2nd observations and used the remaining time in further 
analysis. We also checked the light curve of XIS 1 in the 
energy range of 0.3 to 2.0 keV in the observation periods 
and found no evidence for variation (Fig. 2). 

We constructed instrumental response files (rmfs) and 
effective area files (arf) by running the scripts xisrmfgen 
and xissimarfgen (Ishisaki et al. 2007). To take into 
account the diffuse stray light effects, we used a 20" 
radius flat field as the input emission in calculating the 
arf. We also included in the arf file the degradation of 
low energy efficiency due to the contamination on the 
XIS optical blocking filter. The versions of calibration 
files used here were ae_xil_quantcff_20080504.fits, 
ae_xil_rmfparam_20080901.fits, 

ae_xil_makcpL20080825.fits and 
ae_xil_contami_20071224.fits. Wc estimated the non-X- 
ray-background from the night Earth database using the 
method described in Tawa et al. (2008). 

Wc grouped the spectra to have a minimum number of 
counts in each channel > 50 and used energy range of 0.4— 
5.0 keV in our analysis. This range is broad enough for 
constraining the continuum and also covers the H- and Hc- 
like emission lines of N, O, Ne, Mg, and the L transition 
of Fe. The O VII, O VIII, and NelX lines are clearly visible 
in the spectra (Section 3.2). 

3. Spectral Analysis and Results 

We carried out our data analysis with the Xspec soft- 
ware package, adopting the solar abundances as given in 
Anders & Grevesse (1989). (Hereafter, use of italic type 
indicates Xspec models and their parameters.) Errors 
quoted throughout this paper are single parameter errors 
given at the 90 % confidence level, unless specified other- 
wise. Sections 3.1 and 3.2 give a discussion of our separate 
analyses of the absorption and emission data, while sec- 
tions 3.3 and 3.4 give a discussion of the jointly-analyzed 
data under the uniform and exponential disk models. 

3.1. Chandra X-ray Absorption Spectrum 

We first measured the equivalent widths (EWs) of 
the absorption lines of the highly ionized oxygen ions. 
Bccuase measurement of these narrow absorption lines is 
relevant only to the local continuum, we fit the final PKS 
2155-304 spectrum between 0.55 and 0.7 keV as shown in 
figure 3 using a power-law model modified with absorption 
by the neutral ISM (wabs). The column density of neutral 
hydrogen was fixed to 1.47xl0 20 cm" 2 , which is the value 
determined by the LAB Survey of Galactic HI in this di- 
rection (Kalberla et al. 2005). Three Gaussian functions 
were used to model the OVII K Q , OVIII K Q , and OVII 
K/3 absorption lines (model Al). The measured EWs were 
found to be consistent with those reported by Williams et 
al. 2007. The results are summarized in Table 2 



Once the equivalent widths were determined, we applied 
an absorption line model, absem, to replace the gaussian 
functions in order to probe the properties of the absorbing 
gas. Assuming the temperature and density distributions 
of the hot plasma, the absem model, which is a revision 
of the absline model of Yao & Wang (2005), can be used 
to jointly fit the emission and absorption spectra. (See 
Yao & Wang (2007) and Yao et al. (2009) for a detailed 
description.) For a gas with a uniform density and a sin- 
gle temperature, the diagnostic procedure is summarized 
as follows: (1) A joint analysis of OVII K a and OVII 
Kp directly constrains the OVII column density and the 
Dopplcr dispersion velocity («b)- With the constrained 
Vb, adding the OVIII K a line in the analysis also yields 
the column density of O VIII (model A2). (2) Because the 
column density ratio of O VII and O VIII is sensitive to the 
gas temperature, a joint analysis of the O VII and O VIII 
lines will naturally constrain the gas temperature (model 
A3). (3) Assuming the solar abundance for oxygen and 
given the constrained gas temperature, the O VII (or O VIII 
) column density can be converted to the corresponding 
hot phase hydrogen column density (model A4). Table 3 
gives the results of our fits. The constrained O VII column 
density, 5.9 (+1.2, —0.9) xlO 15 enr 2 is comparable to 
typical values ~ 10 16 cm" 2 obtained from AGN observa- 
tions given in two systematic studies (Fang et al. (2006) 
and Brcgman & Lloyd-Davics (2007)). 
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channel energy (keV) 

Fig. 3. Chandra spectrum of PKS 2155-304 between 0.55 
and 0.7 keV. Fitted model is A4. 

3.2. Suzaku X-ray Emission Spectra 

The Suzaku data were modeled in order to constrain the 
emission measure and the temperature of the halo. For 
this purpose, we first modeled the SXB using a multiple 
component model, since the SXB emission is a superposi- 
tion of such components. We detail this model below. 
3.2.1. Foreground and Background Emission 

We assumed that the SXB consists of four dominant 
components: (1) the Local Hot Bubble (LHB), (2) Solar 
Wind Charge eXchangc in the hcliospherc (SWCX), (3) 
a hot gaseous Galactic halo, (4) the cosmic X-ray back- 
ground emission (CXB; mainly from unresolved extra- 
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Table 2. Spectral fitting results of absorption data with model Al 



Model 




OVII K Q 


VIII K Q 


VII K f3 


Al 


Centroid (eV) 
Sigma (eV) 

Equivalent Width (eV) 


573.81^ 

n 09 + O.25 

U - OZ -0.32 

u.oo+_ 071 


653.1^1 

1 01 +0 - 62 

1 - ui -0.54 

n 077+0116 
U - J ' ' -0.102 


665.8^4 

oi +0 - 97 

u - Ui -o.oi 

n nq +0.058 
U.11J_ 058 



Model Al:wabs(power-law+3x Gaussian) 



Table 3. Spectral fitting results of absorption data with model A2-A4 



Model 


Vb 

(km s- 1 ) 


log 
Abvn 


[Column Density] 
(cm -2 ) 


logT 
(K) 


X 2 /dof 


A2 
A3 
A4 


375+fj 
2901^ 


1 5 7fi+ (J - l)v 

10 - '°-0.08 
1 c: 77+O.O8 
i0 -' '-0.07 


15.56™ 

19.08™ 


R 97+0.02 

°- z '-0.03 
a 90+O.O2 
b - 2 °-0.02 


489.82/474 
498.01/474 
489.84/474 



Model A2,A3,A4: wabs(power-law) X absemx absemx absent 



galactic sources such as AGNs). Because the contribution 
from unresolved Galactic sources is expected to be neg- > 
ligible at high galactic latitudes (\b\ > 30°), we did not =| 
consider such a contribution. The CXB spectrum is well J« 
described by a power-law. c 

In a study of 14 Suzaku blank sky observations, Yoshino g 
et al. (2009) found that there are at least 2 LU (photons "S 
s _1 cm _2 str~ 1 ) of OVII line emission even in those direc- -a 
tions where the attenuation length for the line is less than c 
300 pc. This emission is considered to come from the ^ 
SWCX and LHB, though these contributions are difficult 
to separate with the current CCD energy resolution. After 
Smith ct al. (2007) and Henley et al. (2007), Yoshino et al. 
(2009) found that it can be well represented by a model ^ 
consisting of unabsorbed, optically thin thermal emission ai 
from a collisionally-ionized plasma. The best-fit temper- £ 
ature of this model is log T = 6.06. We therefore use a g 
~ 10 6 K plasma of 2 LU OVII surface brightness as the -S 
SWCX+LHB component. The uncertainty of this esti- ^ 
mate is discussed in section 4.1. a 

Except for the SWCX+LHB component, the observed 
emission has been absorbed by the foreground ISM. In 
the following analysis, we also fix the neutral hydrogen 
column density to be 1.47 x 10 20 cm~ 2 (Kalberla et al. 
2005). 

3.2.2. Spectral Fitting 

To probe the halo gas properties, we used the following 
model to fit our spectra (model El): 
wabs(power - law C XB + vmekal ha io) + mekal LHB +swcx , 
with the photon index of the CXB fixed at 1.4 and with 
the normalization as a free parameter. The temperature 
and the corresponding emission measure (and thus the 
normalization) of the mekalLHB+swcx component were 
set to 1.2 x 10 6 K and 0.0043 pc cm -6 , respectively, 
correspondind to 2 LU of O VII K a line emission. In 
the halo component, we fixed the abundance ratio of 
oxygen to hydrogen to the solar value, and allowed the 
abundances of nitrogen, neon, and iron vary. 

This model fit the spectra from both pointings con- 



- 1 1 1 1 r 




Total 
CXB 

Halo 

SWCX+LHB 
Line Emission - 



rf wt+'t' 




mm 




0.5 1 2 

channel energy (keV) 

Fig. 4. Suzaku spectra between 0.4 and 5.0 keV of Szl 
(top) and Sz2 (bottom) are plotted. Fitted model is E2 
(wabs(power-law CX B + vmekal halo ) + vmekal LH b+SWCX 
+ 3xgaussians). The O and Ne abundance of the vmekalhaio 
(green, dash-dotted) and vmekal^uB+SWCX (blue, dotted) 
are set to be zero and three gaussians (magenta, solid) repre- 
sent Ovn K a , (O vn Kp + O vin K Q ) and Neix K Q 
lines. 



emission 



sistently, except for an apparently higher neon and iron 
abundance in Szl (Table 4) which would be caused by a 
lower temperature of the Szl halo component. It is impor- 
tant to clarify whether this is caused by statistical effects 
or by a true difference in plasma temperature. The sur- 
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Table 4. Spectral fitting results for emission data with the model Ei 



Model Data 


CXB 


LHB- 


-SWCX 




halo 




X'Vdof 




Norm a 


logT 
(K) 


Norm 6 


logT Norm 6 
(K) 


N/O Ne/O 


Fe/O 





El Szl 8.40™ 6.06(fixed) 4.3(fixed) 6.26™ 3.3+^ 6.0+^ 6.5+^ 7.4+f 8 s 148.61/136 

El Sz2 6.45i° : ^ 6.06(fixed) 4.3(fixed) 6.35t° : ^ 3.2^ 4.7tf° l.O^J 147.33/141 

El Szl+Sz2(Szl) 8.30t^ : ^ 6.06(fixed) 4.3(fixed) 6.33^;^ 3.0±^ : | 5.8t^ 3.3±q;1 1.7±J'f 306.82/282 



Szl+Sz2(Sz2) 6.50jlg;lg f t t t t t t 



El t 


Szl- 
Szl- 


hSz2(Szl) 
hSz2(Sz2) 


Q QQ+U.clb 

8 - J °-0.36 


6.06(fixed) 
t 


O.O(fixed) 
t 


6 2^+ U U3 

°- zo -0.02 
t 


4 q+o.v 

*- y -0.6 
t 




t 


4 

t 


313.48/282 


El * 


Szl- 
Szl- 


t-Sz2(Szl) 
hSz2(Sz2) 


fi 47+0-36 


6.06(fixed) 
t 


7.5(fixed) 
t 


6 37+ u - uy 

0,0 ' -0.03 
t 


9 o+0.3 
z ' J -0.3 

t 


t 


o 9+1.3 

°- z -i.o 

t 


1 r+U.9 

t 


299.45/282 



i~ indicates linked parameters 

Szl+Sz2: simultaneous fitting of the data Szi and Sz2 

Model E 1 : wabs(power-lawQx B + nmekal^aio) + mekal^HB+SWCX 



Emission measure of mekal lhb+SWCX is fixed to 0.0043cm -6 which corresponds to 2.0 LU of O VII K Q emission 
T Emission measure of mekalLHB+SWCX is set to as the lower limit 

* Emission measure of mekalLHB+SWCX is set to the upper limit which corresponds to 3.5 LU of Ovu K a emission 

a in unit of photons cm -2 s _1 str — 1 eV — 1 OlkeV 

k Emission measure 10 -3 f n e n v dl: in unit of cm -6 pc 



Table 5. Surface brightness of Ovu, Oviii and Neix 



Model 


Data 


CXB 
Norm a 


Norm h 


halo 
N 


Fe 


OVU K a c 


O VII Kp + 

OVIII K C 


OVIII K Q C 


Ne IX K a c 


xVdof 


E2 
E2 


Szl 
Sz2 


o 9-.+0.62 
°- zl -0.27 
a 07+O.53 
uo ' -0.26 


4 9+0.3 

^• z -0.8 

4 t;+0- 7 

4 -°-0.6 


6.0 (fixed) 
4.7 (fixed) 


7.4 (fixed) 
1.0 (fixed) 


°- uu -0.80 
c 1 c+0.66 
°- io -0.86 


1-45^ 


1 io+ u - 3y 

1 62+ 059 


65+ u - i2 

U.DO_ 26 
U.OO_ 29 


136.39/132 
150.59/137 



model E2: wabs(power-lawc x B + umefeoZ^ a j ) + vmekal^HB+SWCX + 3 X gaussians, where O and Ne abundances of two vmekal are set to 
a in unit of photons cm -2 s _1 str -1 eV — 1 CSfkeV 



Emission Measure 10 3 n e n p dl: in unit of cm 6 pc 
c in unit of LU = photons s _1 cm~ 2 str -1 



face brightness of each line is a better indicator for this 
purpose. 

We next evaluated the surface brightness of O VII and 
O VIII lines by modifying model El (this is model E2). We 
set the O and Ne abundance of the halo and LHB+SWCX 
to zero and used three Gaussian emission lines to represent 
OVII K Q , (O VII K/3 + O VIII Kq,) and NelX K Q emission 
(Fig. 4). Since the XIS resolution is not high enough to 
enable us to distinguish the OVII (656 eV) and OVIII 
K Q (653 eV) lines, they were modeled as a single line. This 
model fitted both spectra with a x 2 /dof of 135.52/132 
and 150.59/137 respectively Assuming the ratio between 
OVII Kg, and OVII K a (= fx) intensities is 0.07 (see foot- 
note 3) , we calculated the O VII, OVIII and NelX surface 
brightnesses as listed in Table 5. Intensities of these lines 
between the two fields are consistent to within the 90% 
confidence level, and we assume that the temperature dif- 

3 H is a slow function of the plasma temperature for thermal emis- 
sion and here the value is 0.056. If the emission is due to SWCX, 
ji = 0.083 (Kharchenko et al. 2003). We averaged these two val- 
ues and used fj, = 0.07 here. See Yoshino et al. (2009) section 
3.1 for details. 



ference is not essential. We plotted the O VII and O VIII 
surface brightness over the Yoshino et al. (2009) results 
(Fig. 5, with 1 a error) and found that the OVII and 
OVIII surface brightness of the PKS 2155-304 direction 
matches the trend of the other 14 fields. 

We next fitted both data sets simultaneously with 
model El by linking parameters of the halo component in 
both observations. The results are shown in Table 5. (Fig. 
6). The emission measure for the model is 3.0 (+0.3, —0.3) 
xlO -3 cm~ 6 pc and the temperature is 2.1 (+0.1, —0.1) 
xlO 6 K. McCammon et al. (2002) reported the emission 
measure and temperature of the absorbed thermal com- 
ponent (=halo) as 3.7 xlO" 3 cm" 6 pc and 2.6 xlO 6 K 
which are comparable to our values. 

3.3. Combined Analysis 

Up to now, we have analyzed the absorption and emis- 
sion data separately and confirmed that the models in- 
cluding the halo component fit both data with a tempera- 
ture of 1.91(+0.09, -0.09) xlO 6 K for the absorption and 
2.14 (+0.15, -0.14) xlO 6 K for the emission spectra. 
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Oyjj Surface Brightness (LU) 

Fig. 5. Relation between Ovn and Ovin surface bright- 
nesses for the 14 (Yoshino et al. 2009) +2 (this work) sky fields 
observed with Suzaku. The horizontal and vertical bars of 
data points show the 1 a errors of the estimate. The contri- 
bution of Ovn emission is corrected for Ovin K Q . The 
diagonal lines show the relation between O vnand O VIII, as- 
suming an offset O VII emission of 2.1 LU and emission from 
a hot plasma of the temperature and the absorption column 
density are shown. The Galactic absorption column density 
of the observation fields are indicated by the maker size of the 
data points. 




channel energy (keV) 

Fig. 6. Suzaku spectra between 0.4 and 2.0 keV. Szl (top) 
and Sz2 (bottom) observations are plotted. Fitted model is 
El (wabs(power-law+vmekal ^aio) J r me ^ a ^LHB+SWCx) an d 
parameters of the halo components are linked in both spectra. 

Assuming that both plasmas are common and uniform, 
the plasma length and density can be calculated using 



the emission measure and the column density. The length 
and density are found to be 4.0 (+1.9, —1-4) kpc and 7.7 
(+2.3, —1-7) xl0~ 4 cm" 3 , respectively. The errors of the 
calculated values are overestimated, since these errors are 
not independent. Moreover, important plasma parame- 
ters such as temperature and velocity dispersion were not 
considered in this simple calculation. 

In this section, using the combined analysis, we will try 
to determine the physical conditions of the halo plasma, 
including the density, the temperature and their distribu- 
tion. 

3.3.1. Uniform Disk Model 

The first step in our combined analysis was to try the 
simplest model: an isothermal plasma with uniform den- 
sity extending up to h kpc above the disk (model CI). 

To perform this combined analysis the emission measure 
and column density have to be linked with a common pa- 
rameter. We chose the equivalent hydrogen column den- 
sity (A^HHot) an d scale height (h) as the control parameters 
and calculated the emission measure. The relation of the 
density n, scale height h, column density -/VH Hot an d galac- 
tic latitude b is described as N-R aat = nh/smb. Thus we 
can use the A4 model for absorption data directly, and re- 
vise the El model to use the vabmkl instead of the mekal 
model. The vabmkl model, an extension of the mekal 
model, was constructed for the combined fit and used the 
column density and plasma length as the fit parameters, 
(see Yao et al. 2009 for a detailed model description) . For 
the halo components of the emission spectra, we fixed the 
abundance ratio of oxygen to hydrogen to the solar value 
and allowed the abundances of nitrogen, neon, and iron to 
vary again. All parameters except for the normalization 
of the CXB components are linked over the two sets of 
emission data. We put lower and upper limits (70-440 km 
s" 1 ) to the velocity dispersion (vb) which represent the 90 
% error range of the values obtained by the absorption 
analysis. 

The model CI fits both data sets (y 2 /dof=802. 78/754) 
and the results are given in Table 6. The column density 
and temperature are consistent with the A4 model (Table 
3), while the temperature and abundance of Ne and Fe 
are not consistent with the El model (Table 4). This is 
because the temperature is mostly constrained by the ab- 
sorption data and the lower temperature for the emission 
spectra preferred the higher abundance to describe the Ne 
and Fe lines. The plasma length is 4.2 (+1.5, —1.2) kpc 
and suggests that under the isothermal assumption the 
halo expands beyond the Galactic disk (~1 kpc). 

3.3.2. Exponential Disk Model 

Observations of edge-on galaxies(ex. Wang et al. 2003, 
Li et al. 2008, Yamasaki et al. 2009) have revealed that 
the intensities of X-ray emission from extended hot gas 
decreases exponentially as a function of height from the 
galactic plane. As a next step in our analysis, we employed 
another simple model to fit the data: an exponential dis- 
tribution model (Yao et al. 2009). In this model, the den- 
sity n and temperature T of the hot gas are distributed 
according to the following equation, 
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Table 6. Combined spectral fitting results with the uniform disk model 



Model 


Data 


CXB 








halo 








X ^/dof 






Norm a 


(cm ) 


h 

(kpc) 


io g r 

(K) 


v b * 
(km s _1 ) 


N/O 


Ne/O 


Fc/O 




CI 


EmissiomSzl 
Emission:Sz2 
Absorption 


a oo+U.39 
8-^-0.38 

D -°'-0.38 


19.08+^ 

t 
t 


4.2+^ 
t 


D - Z '-0.02 
t 
t 




t 


— 1.5 

t 


t 


802.78/754 



i~ indicates linked parameters 

model CI: wabs(power-law+vmekal)+mekal for the emission, wabs(power)x(absem)' i for the absorption 
*Parameter range is limited to 70-440 km s — 1 
a in unit of photons cm" - 2 s- 1 sir' 1 cV- 1 OlkeV 



n = n e- z / h ^, T = T Q e^ z ^ hT ^ 1 1 = h T /h n (1) 

where Z is the vertical distance from the Galactic plane, 
no and To are the density and temperature at the plane, 
and h n and hx are the scale heights of the density and 
temperature, respectively, and £ is the filling factor, which 
is assumed to be 1 in this paper. Thus the equivalent 
hydrogen column density of the hot gas (-/VH Ho t) i s ca l~ 
culated as Nn Hot = J Q ndl = f n cxp(—Z/h n )dZ/sinb = 

The models vabmkl and absem can also be used in an 
exponential disk model using the additional parameter 7 
(see Yao et al. (2009) for detailed description). We there- 
fore used the same model as used in the uniform model 
here (model C2). For fit parameters, for convenience we 
used the column density -/Vhhoi instead of no- 

We jointly fitted the emission and absorption data using 
this exponential disk model. The parameters obtained are 
summarized in Table 7. We first fixed the velocity disper- 
sion (vb) at 290 km s — 1 . We next examined the robust- 
ness of the temperature (To), column density (-/Vh hot ), 
and scale height (h n ), as a function of 7, Vb, and the in- 
tensity of foreground SWCX intensity. We found that all 
parameters are consistent to within 90% statistical errors. 
When we fitted with Vb allowed to vary freely, the best-fit 
value of vi, became km s _1 . Though this is above 

the thermal velocity (~ 30 km s _1 ), it is a smaller value 
than that obtained from the absorption spectrum which 
determined the ratio between the OVII K Q and Kp lines. 
In the exponential disk model, low (3 x 10 5 K < T < 10 6 K) 
temperature plasma can exist in the outer regions, which 
contribute only to the OVII absorption line. This might 
cause the smaller Vb value. The cooling time of such low 
temperature plasmas is very short, and the actual situ- 
ation will not follow such a simple exponential model in 
this temperature range. We therefore fixed Vb at 290 km 
s , as the best-fit value from the absorption analysis. 

Confidence contours of h n , To and A^j Hot versus gamma 
are plotted in figure 7, over-laid on those of the LMC 
X-3 direction (Yao et al. 2009). We then obtained the 
scale height for the temperature gradient as h t = 5. 61^2 
kpc and the gas density at the galactic plane as n n = 
(1.4+ o\l) x 10~ 3 cm" 3 (Figure 8). This values is typi- 
cal for the mid-plane plasma density (Cox 2005). As the 
high temperature plasma close to the Galactic plane can 
emit Fe and Ne lines efficiently, the spectrum can be fitted 
without an abundance of heavy element higher than the 




Fig. 7. 68%, 90%, and 99% confidence contours of h n , T , 
and A r jj Hot vs. 7, obtained from the combined fits to the 
X-ray absorption and emission data. Colored thick lines are 
for the PKS 2155-304 sight line, while the black thin lines are 
for the LMC X-3 sight lines (Yao et al. 2009). In the panel 
(a) the scale height of the temperature (/it) is constant along 
the dashed lines. 

solar value. The emission weighted temperature calcu- 
lated with best fitted parameters using the intensity ratio 
of O VIII to OVII becomes 2.2(+0.1, -0.1) xlO 6 K. 
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Table 7. Combined spectral fitting results with the exponential disk model 



Model 


Data 


CXB 








halo 










X^/dof 






Norm a 


(cm - ) 




10 g ^0 




7 


NT /O 


Mr. /fl 

in e/ \j 


re/ w 










(KpC) 




(km s ) 












C2 


TV, ICJGir^Tl " ^ r 7] 
JJjlIllSSlUll. U£il 


8 - ZD -0.37 


iy.iu_ 07 


n q+u.y 

Z - J -0.8 


D - 4U -0.05 




9 44+1.11 


r o+l.B 
°-°-1.3 


o 1+1.6 
6 - L -l .2 


i 5 +1 - u 

1 '°-0.7 






Emission: Sz2 




t 


t 


t 




t 


t 


t 


t 






Absorption 




t 




t 


290 (fixed) 


t 








792.76/757 


C2 


EmissiomSzl 


8.20+"; 39 


19 13+ 007 


9 9 +0.5 
z - -0.7 


fi 48+0.04 




l.O(fixed) 


6 1+ 1 - 8 

D.l_ 1 4 


2 4 +a9 
z '^-0.9 


i o +0 - 6 






Emission: Sz2 


6 40+ 38 


t 


t 


t 




t 


t 


t 


t 






Absorption 




t 




t 


290 (fixed) 


t 








795.64/758 


C2 


Emission:Szl 


8 25+ 033 
°- zo -0.38 


19 10+ 007 

±M.±U_ 07 


2 4+ - 9 


f- OO+0.02 




3.5(fixed) 




o q+1.2 
°-°-0.8 


-, 7 +0.3 
1 - '-0.5 






Emission:Sz2 


D -^°-0.37 


t 


t 


t 




t 


t 


t 


t 






Absorption 




t 




t 


290 (fixed) 


t 








793.64/758 



C2 


Emission:Szl 
Emission:Sz2 
Absorption 


a 17+0.37 

6 37+ 37 
-0.38 


19 41+ 019 

t 
t 


t 


6 ^1+ 016 
b - si -o.io 

t 

t 


70 (fixed) 


n 43+ 116 
u '*°-0.23 

t 

t 


= 7+1-5 
t 


2 3+ 1 - 

t 


1.0+" 

— 0.5 

t 


789.75/757 


C2 


Emission:Szl 
Emission:Sz2 
Absorption 


o 97 + 0.38 

°"''-0.33 
6 47+0.38 
D '^'-0.32 


lg -,0+0.06 
13 - 1J -0 .07 

t 

t 


I q+0.5 

t 


6 40+ - 04 

°-^ u -0.05 
t 
t 


440 (fixed) 


2.84ti; 3 4 

t 
t 


t 


3 2+ 1 ' 1 
°- -1.0 

t 


1 6 +0 - 5 
^•"-o.s 

t 


795.44/757 


C2t 


Emission:Szl 
Emission:Sz2 
Absorption 


o O1+0.33 
8->5l_o.4l 

6 ^2+ 33 


19.08i° ;« 

t 
t 


i c+0.6 
i -°-0.5 

t 


6 36+ 03 

t 
t 


290 (fixed) 


°'° -1.97 

t 
t 


4 9+1-1 
t 


o 9+0.7 

t 


2 0+°- 9 
^• U -0.8 

t 


803.71/757 


C2* 


EmissiomSzl 
Emission:Sz2 
Absorption 


8 lfi+0.38 
°- ib -0.39 
6 3K+0 38 
°- m -0.3S 


19 16+ ' 09 

1».1D_ og 

t 
t 


t 


6 51+ ai0 
D - o± -0.07 

t 

t 


290 (fixed) 


1 47+0.52 

t 
t 


7 9+2.4 
'• z -1.8 

t 


9 9 + 1.3 
z - -1.0 

t 


9+°- 6 
t 


783.37/757 



f indicates linked parameters 

model C2: wabs(power-law+vabmkl)+mekal fot the emission, wabs(power)x(absem) 3 for the absorption 
tEmission measure of mefcai£//a_|_5vy C 7x is se t to as the lower limit. 

^Emission measure of mekal^ug^gyvcx ls se t to upper limit which corresponds to 3.5 LU Ovn K a emission 
a in unit of photons cm -2 s~ 1 str -1 eV -1 @lkeV 



4. Discussion 

4-.1. Uncertainty due to of LHB and SWCX 

Because our knowledge about the temporal and spatial 
variability of the SWCX and the LHB is limited, there 
are uncertainties due to the assumption of their intensity. 
These uncertainties could result in large uncertainties in 
our results. 

To assess this uncertainty, we estimated the lower and 
upper values of the LHB and SWCX contributions and 
evaluated the parameters of the halo components again. 
The lower limit of the contribution is zero. As for the 
upper limit, we adopt 3.5 LU for the OVII emission, as 
obtained by the MBM 12 shadowing observation (Smith 
et al. 2007). As the heliospheric SWCX is caused by the 
collision between the Solar wind and the neutral ISM, the 
estimated emissivity has a peak around the ecliptic plane 
(Koutroumpa et al. 2007 , Lallcmcnt et al. 2004). MBM 
12 is located at (A,/3)=(47.4, 2.6) in ecliptic coordinates, 
while PKS2155-304 is at (A,j8)=(321.2, -16.8). Thus we 
assume that the heliospheric SWCX contribution in the 
PKS 2155-304 direction could not be larger than that for 
MBM12. 

The results using these lower and upper limits are shown 
in Table 4 and Table 7. Though the best fit values are 



slightly changed, they are consistent with the previous 
analysis. 

4-2. Comparison with the Results for LMC X-3 

We compared our results with those of the LMC X-3 
direction, as is summarized in Table 8. The directions of 
the LMC X-3 and PKS 2155-304 are (l,b) = (273.6,-32.1) 
and (17.7,-52.2). The fact that we obtained similar val- 
ues for the two directions indicates that the hot halo is 
common in the big picture and can be explained with the 
exponential model of the column density, scale height and 
temperature as ~ 2 x 10 19 cm -2 , a few kpc and ~2x 10 6 
K. As the distances to the targets are 50 kpc for LMC X- 
3 and 480 Mpc for PKS 2155-304, the consistency of the 
parameters of the exponential disk suggests that there is 
little contribution from beyond LMC X-3, or from a very 
extended halo of a 100 kpc scale. 

4-3. Distribution of the VII and O VIII Emitting/ 
Absorbing Gas and Its Origin 

We calculated the distribution of O VII and O VIII ions 
and their emissivities assuming the best fit parameters at 
7 = 2.44 and at 7 = 1.0 and 3.5 (Fig. 9). 

We then estimated the total radiative energy loss from 
the thick disk distributed exponentially. Assuming solar 
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Table 8. Disk model parameters for two sight lines 



Direction 


!0g ^H H ot 

(cm - ") 


(kpc) 


logT Q 
(K) 


7 


Ne 


Fc 


PKS 2155-304 
LMC X-3t 


19.10±H 8 r 
19.36±8;|f 


O+0.9 
z -°-0.8 

2.8t 3 ; 6 8 


fi 40 +0 - uy 
fi 5fi+ 011 

°- OD -0.10 


2.44±iii 


1+1.8 

1 7+O.6 
' -0.4 


1 5 +1 ' U 

1 '°-0.7 

o q +0 - 2 

u.y_ 2 



f from Yao et al. 2009 
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Fig. 8. 687o, 90%, and 997c confidence contours of T and 
iVjj Hot vs. scale height h n obtained in the joint fit to the 
X-ray absorption and emission data. In the upper panel the 
density at the plane no is constant along the solid and dashed 
lines. 

abundances, best fit parameters and ionization fraction 
and emissivity as taken from SPEX 4 , we obtained the 
energy loss rate as a function of the distance from the 
Galactic plane Z (Fig. 10). We then integrated the en- 
ergy loss rate until the temperature of the exponential disk 
become lower than 10 5 5 K. Because our results are based 
on X-ray observations, it is difficult to detect plasma of T 
< 10 5 ' 5 K. We obtained a total radiative energy loss rate of 
7.2 x 10 36 erg s~ x kpc" 2 in 0.001-40 keV and 1.8 x 10 35 erg 
s" 1 kpc -2 in 0.3-8.0 keV. These values are consistent with 
the X-ray luminosity of other spiral galaxies (Strickland 
et al. 2004). 

We next compared the energy loss rate with the energy 
input rate from SNe. According to Ferriere (1998), the SN 

4 http:/ /www.sron.nl/indcx.php?option=com_contcnt 
&task=view&id=125&Itemid=279 



Fig. 9. The density of Ovu and Ovm ion(top) and the 
emissivity of Ovu and Ovm lines (bottom) as a function 
of the height from the galactic plane under the best fit pa- 
rameter of 7=2.44 (solid line), 7=1.0 (dashed line) and 7=3.5 
(dash-dotted line). 

rate near the sun is 19 Myr" 1 kpc~ 2 for type II SNe and 
2.6 Myr -1 kpc -2 for type la SNe, respectively. Assuming 
each SN explosion releases 1 xlO 51 ergs, the total input 
energy is then 7 xlO 38 ergs s" 1 kpc" 2 . If 1 % of the SN 
explosion energy is input to the hot halo, the total energy 
loss can be compensated. 

4-4- Consistency with O VI Absorbing Gas 

It is not clear that our model is consistent beyond ~ 5 
kpc where the temperature of the gas is below ~ 10 6 K 
and O VI ion becomes dominant. 

Williams et al. (2007) found two local Ovi absorption 
lines in the FUSE PKS 2155-304 spectrum and reported 
column densities of 1.10±0.1 x 10 14 and 8.7±0.4 x 10 13 
cm" 2 . Our exponential disk model expects OVI column 
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Fig. 10. Radiative energy loss rate (red, solid) and cool- 
ing time (blue, dashed) as a function of the distance from 
the Galactic plane. The temperature is indicated by the 
solid black line. The emissivity is calculated from the mekal 
model, using a script made by Sutherland, 
(http: / /proteus. pha.jhu.edu/~dks/Code/ 
Coolcurve_create / index, html) 

densities of 3.8 x 10 13 , 1.4 x 10 14 , and 2.1 x 10 13 cm" 2 
with the best fit parameters when 7=2.44, 1.0, and 3.5 
respectively. 

Howevera plasma emitting O VI lines cools very rapidly 
and it would be difficult to maintain such plasma existing 
high above the Galactic plane. Radiative cooling is accel- 
erated by the density fluctuations. Thus OVI absorbing 
gas can be a patchy or blob-like condensation. To discuss 
this problem, energy and matter flow models are needed, 
which is beyond the focus of this paper. 

5. Summary 

We have analyzed high resolution X-ray absorp- 
tion/emission data observed by Chandra and Suzaku to 
determine the physical state of the global hot gas along 
the PKS 2155-304 direction. 

1. Suzaku clearly detected OVII K a , OVIII K Q and 
OVII Kp lines. The surface brightnesses of OVII 
and O VIII in this direction can be understood in the 
same scheme as obtained by other 14 observations 
(Yoshino et al. 2009). 

2. By the absorption analysis, column density is mea- 
sured as 3.9 (+0.6,-0.6) cm -3 pc and temperature 
is measured as 1.91 (+0.09,-0.09) xlO 6 K. By the 
emission analysis, emission measure is measured as 
3.0 (+0.3,-0.3) xl0~ 3 cm~ 6 pc and temperature is 
measured as 2.14 (+0.15,-0.14) xlO 6 K. 

3. Combined analysis using the exponential disk model 
gives a good fit with x 2 / do1 01 789.65/756 to 
both emission and absorption spectra. The gas 
temperature and density at the Galactic plane 
are determined to be 2.5(+0.6, —0.3) x 10 6 K and 
1.4(+0.5, -0.4) x 10 -3 cm -3 and the scale heights 
of the gas temperature and density 5.6(+7.4, — 4.2) 
kpc and 2.3(+0.9, — 0.8) kpc, respectively. 



4. The results obtained by the combined analysis are 
consistent with those for the LMC X-3 direction. 
This suggest that the global hot gas surrounding 
our Galaxy has common structure. 
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